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REMOTE SENSING OF SOIL MOISTURE 
WITH MICROWAVE RADIOMETERS-II 


T. Schmugge, T. Wilheit, W. Webster, Jr., and P. Gloersen 

Goddard Space Flight Center 
Greenbelt, Maryland 


INTRODUCTION 

Interest in developing techniques for the remote sensing of soil moisture has increased 
significantly in recent years in response to the need for this type of data covering large 
areas. The information is needed in the fields of meteorology, hydrology, and agricul- 
ture. In the field of meteorology, knowledge of soil moisture is necessary for the large- 
scale assessment of the moisture flux into the atmosphere. Knowledge of the moisture 
content in a surface layer about 10 cm thick is considered important for predicting the 
runoff following a rain, and the successful germination and early growth of plants. However, 
to predict crop yields, the moisture content in the root zone must be known. This depth 
is approximately 1 m and is probably not directly measurable by remote sensing. In addi- 
tion, remotely sensed ! soil-moisture data would be of use in locating potential outbreaks 
of certain insect pests which deposit their eggs in the soil, since the subsequent develop- 
ment of the eggs is dependent on soil-moisture content (Reference 1). 

Recent efforts have been made to correlate soil moisture with several parameters that can 
be remotely sensed, such as surface albedo (Reference 2), the diurnal range of surface 
temperature (Reference 3), microwave backscatter coefficient (Reference 4), and micro- 
wave emissivity (Reference 5). 

The variation of the microwave emissivity for a soil has been studied extensively by truck- 
mounted radiometers in fields (References 6 and 7). These measurements have indicated 
the sensitivity of microwave emissivity to soil moisture and the dependence of this sensitivity 
on surface roughness and vegetative cover. Researchers have concluded that the longer 
wavelength radiometers, for example 21 cm, are preferable for the remote sensing of soil 
moisture. Radiometers operating at this wavelength showed lower sensitivity to the effect 
of surface roughness and a greater capability to sense moisture variations through a vegeta- 
tive canopy than those operating at shorter wavelengths. 

Aircraft observations of these effects have been made both in this country by Schmugge 
et al. (Reference 5) and in the Soviet Union by Basharinov et al. (Reference 8). The 
research performed in the United States determined a nonhnear dependence of microwave 



brightness temperature (Tg ) on soil moisture from aircraft platform observations. Little 
change in brightness temperature was observed at a wavelength of 1 .55 cm for soil-moisture 
values less than 10 or 15 percent (depending on soil type), while above these values, a linear 
decrease of about 3 K/1 percent soil moisture was observed. The soil-moisture values re- 
ported were the average moisture contents in the top 15 cm of the soil, and 'these results in- 
dicate that the nonlinear behavior was due to the drying out of a thinner surface layer. 
Calculations using a radiative transfer model for layered dielectrics indicated that the ob- 
served brightness temperatures were determined by the moisture content of a surface layer 
only ~0.1 of a wavelength thick. In this paper, results are given from flights with the NASA 
Convair-990 aircraft in 1972 and 1973 which confirm that the emissivity is indeed determined 
by the moisture content of this layer. 

BACKGROUND 

The intensity of the emitted radiation is essentially proportional to the product of the temp- 
erature and emissivity of the surface. This product is commonly referred to as brightness 
temperature (Tg ). The radiation received by an airborne or spaceborne radiometer at an 
altitude (H) above the groimd is given by 

Tg = T (H, 0) [R TgjjY (1 “ R) Tgug^p] + ^ ) 

where r is the atmospheric transmission and R is the surface reflectivity. A schematic is 
presented in figure 1 . The three terms of the equation represent the three sources of the 
observed radiation: first, reflected sky brightness (Tgj^y ); second, the thermal emission from 
the surface (Tg^gp); and third, the emission from the intervening atmosphere (T^^m 
factor r is atmospheric absorption of the radiation coming from the surface. For the low- 
altitude observations reported in this paper, the atmospheric effects are negligible as will be 
shov/n later. Thus, we can concentrate on the surface term. The quantity (1 - R) is the 
emissivity for the surface. The large variation at microwave wavelengths of this quantity for 
terrestrial surfaces causes the major variations in the observed brightness temperature. 

The large dielectric constant of liquid water produces most of this variation in the surface 
emissivity. The real and imaginary parts of the dielectric constant of pure water at 20°C 
are plotted in figure 2 for the 100- to 0.1 -cm wavelength range. These data are based on the 
measurements of Lane and Saxton (Reference 9). In this region of the spectrum, the real 
part decreases from about 80 at 100 cm to less than 10 at 0.1 cm. The imaginary part 
reaches its maximum around 1 .5 cm where the real part is changing most rapidly. 

Most dry soils have dielectric constants of about 3 for the real part and less than 0.0 1 
for the imaginary part. The addition of water causes marked changes in the dielectric 
properties. Figure 3 contains data on the variation of these dielectric properties with 
moisture content at two wavelengths— 21 cm (Reference 10) and 1.55 cm. These measure- 
ments were made on two soils having very similar textures, that is, particle-size distribution. 
Thus, the differences between the two sets of measurements arise from the wavelength 
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difference. T his results from the differences in the dielectric properties of water at these 
two wavelengths, that is, for wet soils the real part is larger at 21 cm while the imaginary 
part is larger at 1 .55 cm. The resultant emissivities have been calculated, assuming a uni- 
form dielectric medium, from the Fresnel equation (Reference 11). The range of emis- 
sivity from dry to wet was 0.94 to 0.50 at 21 cm and 0.95 to 0.60 at 1.55 cm. Thus, 
observations at the longer wavelengths would have greater sensitivity to soil-moisture varia- 
tions. 

It can be further seen in figure 3 that the addition of water has very little effect on the 
dielectric properties of the soil at low-moisture contents (<10 percent). Presumably this 
is due to the interaction of the water molecules with the soil particles which reduces the 
polarizability of the water in a thin layer around each particle. The frequency variation 
of the dielectric constants for a silty clay soil with a similar moisture content has been 
observed by Hoekstra and Delaney (Reference 12). The curves have shapes similar to 
those of water (figure 2), but the peak of the imaginary curve occurs at a|lower frequency, 
between 1 and 2 GHz. They were able to fit the soil data with a modified Debye equation 
containing two terms. One term has a characteristic relaxation frequency of 1 .5 GHz 
and is representative of the bound water, and the second has a relaxation frequency of 
20 GHz which is roughly that of free water. As the water content of the soil increases, 
one would expect that the free-water term would become dominant. Unfortunately, 
Hoekstra and Delaney were not able to perform their measurements at water contents 
greater than 1 5 percent by weight. Thus, they did not get into the region of moisture con- 
tent for clay soils where the soil water has its largest dielectric effect (the region of steepest 
slope in figure 3). 

At 1.4 GHz, Lundien (Reference 10) has observed that the dielectric constant versus the soil- 
moisture curve can be approximated by two straight lines, and that their intersection is a 
linearly increasing function of the fraction of particles finer than 0.02 mm in the soil. The 
intersection for a sand was at 7 percent and for a heavy clay at 21 percent. Thus, the di- 
electric properties of a soil are dominated by the amount of free water in the soil, which 
in turn is related to the soil texture. 

The lines in figure 3 are linear regression fits to the data in the region represented by each 
line. These regression results will be used in calculations discussed later in this document. 

Since the thermally emitted microwave energy originates within the volume of the soil 
and not at the surface, the term (1 - R) Tg^j^p is a simplification of the actual situation be- 
cause it assumes a uniform temperature throughout the emitting medium. Therefore, this 
term should be an integral over the entire thickness of the soil under consideration, that is. 


(1-R)Tsurf = (1-R) 



T(z)a(z) exp 




dz 


( 2 ) 
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where T(z) and a(z) represent the subsurface temperature and absorption coefficient pro- 
files, respectively. WiUieit (Reference 13) has developed a model in which the integral is re- 
placed, by a sum over many layers: 


(1 - R) Tgyj^p fj Tj (3) 

where f. is the fraction of the radiation incidence on the air-soU interface that would be 
absorbed in the ith layer ,and Tj is the temperature of this layer. The values of f. are deter- 
mined by applying the 'electromagnetic boundary conditions to determine the electric fields 
above and below each interface. These fields are then used to determine the energy fluxes 
entering and leaving each layer. The computations indicate that the radiation from the soil 
is characterized by two sampling depths; reflective and thermal. The reflectivity is charac- 
terized by changes in the real part of the index of refraction over a sampling depth: 8^ ~ 0.1 X, 
where X is the wavelength in the medium. The thermal sampling depth is determined by the 
losses deeper in the medium, and is given by 



where x. is the depth of the ith layer. For a uniform dielectric this reduces to 


S 


T 




An Im(n) 


(4) 


(5) 


For a low-loss dry soil, 5.^. will be an order of magnitude larger than 8^., while for a wet soil, 
it will be only sUghtly larger. 


Another parameter of interest is the average soil temperature over this thermal sampling which 
is given by 



( 6 ) 


where T. is the physical temperature of the ith layer. The ratio of Tg to T^^ is an effective 
emissivity for the soil. 


The soil-moisture information of the most interest is the amount of water in the soil which 
is not tightly bound to the soil particles. This water can freely evaporate from the soil and 
is readily absorbed by plants. The data presented in this report indicate that the emissivity 
of soil is lowered only when the soil water is not bound. The reason for this is that the tightly 
bound water does not have the same dielectric properties as free water and, more importantly, 
it can only move into the surface layer of the soil very slowly. This movement is largely in 
the vapor phase and occurs as a result of alternate condensation and evaporation on the 
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particle surfaces. Only when moisture is in the surface layer of the soil can it affect the 
microwave emission from the soil. 

The state of the water in the soil is described by the pressure potential. At low-moisture 
levels, the pressure potential is the tension with which water is held by soil particles. In 
the intermediate range, the pressure potential is determined largely by the radiilof curvature 
of water films between soil particles. In figure 4, representative plots of the relation between 
volumetric water content and pressure potential are presented (Reference 3). The zero po- 
tential level is the saturated soil situation. The locations of the -1 /3-bar and -15-bar pressure 
levels are indicated because they are frequently taken to be the pressure levels for the field 
capacity (FC) and wilting point (WP) conditions of the soil. This convention will be followed 
here. The amount of water in the soil at field capacity is that which remains in a soil two 
or three days after having been saturated and after free drainage has practically ceased. As 
the name implies, the wilt point is the moisture level at which plants experience difficulty 
drawing water from the sod. Thus, the FC and WP soil-moisture values give a quantitative 
measure of the water-holding capacity of a soil. The difference between the two is the 
available water capacity in the sod. As the curves in figure 4 indicate, FC and WP depend 
on soil type. The values for the four soils are presented in table 1 . 


Table 1 

Moisture Content of Selected Sods at 
WP, FC, and Saturated Conditions 


Soil 

WP at -15 bar 
(cm^/cm^) 

FC at -1/3 bar 
(cm^/cm^) 

Saturation at 0 bar 
(cm^/cm®) 

Navajo Clay 

0.22 

0.35 

0.70 

Cashion Silty Qay 

0.22 

0.33 

0.50 

Avondale Loam 

0.11 

0.25 

0.44 

Gran Sandy Loam 

0.06 

0.15 

0.36 


Thus, the avadable water capacity of the clay soils is no greater than that of the Avondale 
Loam, even though the magnitudes are much greater. 

Since it would be very difficult to obtainj curves like those in figure 4 for the sods in aU 
the sampled fields, we attempted to relate FC and WP to the sod textures of the sampled 
fields. This was based on the work of Salter and Williams (Reference 14) who studied the 
use of regression analysis to relate particle-size composition (sod texture) to the available 
water capacity for a sod. They concluded that the moisture characteristics could be cal- 
cidated from texture information with reasonable accuracy; that is, the upper and lower limits 
of avadable water capacity could be estimated to within 1 0 or 20 percent of the measured 
values. Therefore, a multiple linear regression analysis was made on 100 sets of soil tex- 
tures and moisture characteristics, that is, the moisture contents at the -l/3-bar>(FC) and 
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-15-bar (WP) potentials. These measurements were made on soils from the Phoenix area 
(Private communication, Phoenix Soil Conservation Office, 1974) and from the Rio Grande 
Valley of Texas (Reference 15). The range of textures included in the regression covered 
those observed in the fields. The results of the regression for WP, expressed in weight per- 
cent, were 

WP = 7.2 - 0.07 X SAND + 0.24 X CLAY (7) 

where SAND and CLAY represent their respective soil fractions in percent. The multiple 
correlation coefficient for this regression was 0.945. The regression results for FC are 

FC = 25 .1 - 0.21 X SAND + 0.22 X CLAY (8) 

with a multiple correlation coefficient of 0.904. The coefficient of variation (standard 
estimate of error divided by the mean) was 0.15 for both of these regressions. The moisture 
characteristics of a soil depend on many factors in addition to soil texture, such as bulk den- 
sity of the undisturbed soil and percent organic matter, but texture (sand, silt, and clay frac- 
tions) was the only parameter that could easily be determined for all of the soils involved. 

It is presumed basing the regression on the actual field soils used in the soil surveys adequately 
takes these factors into account. 

The effect of a vegetative canopy over the soil will be that of an additional dielectric layer, 
the mixture of plant matter and air, between the radiating soil arid the radiometer. Sibley 
(Reference 16) has estimated the dielectric constant for this layer using a Weiner mixing 
formula and found that a canopy with a height density product of 4 cm (that is, 20-cm high 
wheat with 20 percent of the volume containing plant matter) increased the brightness 
temperature above that expected for a bare, moist soil by 5 K at 1 GHz, 40 K at 3 GHz, and 
70 K at 5 GHz. Field measurements at 1.42 GHz (Reference 7) have qualitatively supported 
this result. 

EXPERIMEIMTAL DETAILS 

The flights described in this report were flown in March 1972 and February 1973. The 
aircraft used in the experiment was the NASA Convair-990. The test sites were the agri- 
cultural areas around Phoenix, Arizona and Imperial Valley, California. The times of the 
individual flights and the meteorological conditions at each site are described in table 2. 

The altitude for all the data flights was 0.6 km above the ground with a ground speed of 
1 10 to 120 m/s. At this altitude, the effects of the atmosphere on the observed brightness 
temperatures (see equation 1 ) are minimal for all of the wavelengths considered in this re- 
port, as indicated in table 3. 

These calculations were performed using an atmospheric profile obtained from the nearest 
rawinsonde station (Tucson). This was a dry atmosphere with 12-percent relative humidity 
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Table 2 

Description of Data Gathering Flights 


Date 

Site 

Local 

Standard Time 

Meteorological Conditions 

March 7, 1972 

Phoenix, Arizona 
100 fields 

14:04 to 14:45 

Clear, 100 percent sunshine 
Air Temperature = 3 1°C 
1 7 percent Relative Humidity 

March 11, 1972 

Phoenix, Arizona 

12:25 to 13-.J2 

0.3 Cirrus layer at 7600 m 
Air Temperature = 31 .4°C 
10 percent Relative Humidity 



15:35 to 16:23 

0.5 Cirrus layer at 7600 m 
Air Temperature = 32°C 
10 percent Relative Humidity 

March 11, 1972 

Imperial Valley 
50 fields 

11:00 to 11:50 

0.5 Cirrus at 7600 m 

Air Temperature = 28°C 

2 1 percent Relative Humidity 



17:09 to 17:55 

0.5 Cirra stratus at 6700 m 
Air Temperature = 29°C 
1 8 percent Relative Humidity 

March 13, 1972 

Phoenix, Arizona 

13:10 to 13:55 

0.5 Cirrus at 7600 m 

Air Temperature = 29°C 

1 1 percent Relative Humidity 

February 8, 1973 

Phoenix, Arizona 
80 fields 

13:45 to 14:35 

0.2 Cumulonimbus at 1700 m 
Air Temperature = 19°C 
42 percent Relative Humidity 

February 8, 1973 

Imperial Valley 
50 fields 

15:28 to 15:47 

No Data Available* 


♦The closest station to the Impeiial Valley, Yuma, Arizona, was not operating in 1973, and therefore data are not 
available. 


at the surface and a total of 1.1 cm of water in the column. The AT values are the magnitude 
of the atmospheric effects for a surface reflectivity of 0.3 which would be the approximate 
value for a wet field. For a lower reflectivity, the AT values would be even smaller. Thus, 
it is apparent that atmospheric effects are small and can be neglected. 
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Table 3 

Atmospheric Effects 


Wavelength 

(X) 

(cm) 

■ 

Sky Brightness 
(K) 

Transmissivity 

(r) 

(K) 

Atmospheric 

Contribution 

^^ATM ^ 

(K) 

Net Effect 
(K) 

21 

•5.0 

0.999 

0.3 

1.6 

6 

5.9 

0.998 

0.5 

1.9 

2.8 

6.6 

0.997 

0.8 

2.1 

1.55 

11.7 

0.991 

2.7 

3.8 

0.8 

19.8 

0.985 

4.5 

7.4 


There were two tracks flown at each site: In Arizona there was a north-south track west 
of the city of Phoenix, and an east-west track just south of the city. In Imperial Valley 
there were two east-west tracks, one at the north end of the Valley over the town of 
Calipatria and one at the south end over the town of Imperial. The latter track was re- 
placed for the 1973 mission with a second one at the north end of the Valley to simpli- 
fy ground truth acquisition. For the 1972 flights, ground data were acquired for approx- 
imately 100 16-hectare (ha) fields around Phoenix and 50 such fields in Imperial Valley. 

The details of the 1972 ground observation program are described elsewhere (Reference 17). 
To ensure that data were obtained over the middle of these fields, two passes were made 
over each flight line and the best data were used in the analysis. A different approach 
was taken for the 1973 flight so that only fields that were 800 m wide and 400 m long 
(area = 32 ha) were sampled, and there was only one pass over each flight line. 

The average soil moisture for a field was determined for depth intervals of 0 to 1 cm (where 
possible), 0 to 2.5 cm, 0 to 5 cm. and 0 to 10 cm (1972 flights only). The averages were ob- 
tained using samples from four locations in each field. The four samples for each depth 
were combined and the resultant soil moisture was determined gravimetrically. The soil- 
moisture values presented in this document are expressed on a dry-weight basis. The sampling 
in Imperial Valley was done by the John D. Hess Testing Corporation, and around Phoenix, 
by personnel from the Salt River Project. For the 1973 program, a textural analysis was 
performed on the soil sample obtained from each field. The soils from Imperial Valley were 
found to be predominantly clay loams and clays with a few sandy loam fields at the eastern 
edge of the Valley, while the soils from the Salt River Valley were generally not as heavy and 
were more evenly distributed within the range between sandy loam and clay loam. This 
knowledge of the soil texture enabled us to make estimates of the moisture characteristics 
for the soils using the regression results discussed previously in equations 7 and 8. 

For the 1972 flints, textures of the soil samples were not determined but were estimated 
by the field crews. In those fields sampled in both missions, the agreement between the 
laboratory determinations and the field estimates was very good. 
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The characterization of the soil moisture for a 16- or 32-ha field by the average of only 
four samples introduces considerable uncertainty. The problem of the variability in soil- 
moisture estimates has been studied by Reynolds (Reference 18), who found that, “A sample 
size of 10 is usually sufficient to estimate the moisture content with a 0.1 coefficient of var- 
iation at the 95-percent probability level.” While these measurements were made on 3.1 X 1.9-m 
plots for 0- to 8-cm depths, Reynolds states that the sample sizes can be taken as representa- 
tive of much larger areas. More pertinent to our situation is an experiment performed by 
Schmugge et al. (Reference 19) on a 16-ha field south of Phoenix in which an estimate of the 
accuracy of the four-point measurement for a field of this size was made. This estimate is 
based on 36 uniformly spaced samples which were divided into four quarters with nine loca- 
tions. The mean for a sample consisting of a data point from each quater of the field was cal- 
culated 100 times. The distribution of results was then compared with the overall field mean. 
The results for the limits of conficence at the 90-percent level are listed in table 4. 


Table 4 

Limits of Confidence on Ground Measurements 


Depth 

(cm) 

Field Mean 
(percent) 

Limits 

(percent) 

Oto 1 

15.8 

-2.9 to +3.5 

1 to 2 

18.4 

-2.6 to +3.1 

2 to 5 

21.8 

-2.5 to +2.4 

5 to 9 

23.1 

-1.8 to +2.6 

9 to 15 

23.9 

-1.5 to +1.7 


Thus, for the 0-to 1-cm layer sampled in this experiment we would expect that our measured 
values are within 20 percent of the actual value. 

Another source of uncertainty is the time separations between the aircraft overflight and the 
ground measurements. Jackson et al. (Reference 20) have observed that there is a large di- 
urnal variation in the moisture content of the surface layers, for example, ~1 5 percent for 
the 0- to 1-cm layer three days after irrigation. The aircraft overflights occurred between 
noon and 3 p.m., while the ground acquisition took place over a longer time interval— all day 
at Imperial Valley. Therefore, those measurements made early in the morning may be 5-to 
10-percent wetter for the surface centimeter than was actually the case at the time of over- 
flight for a wet field. Fields sampled in the afternoon should be in good agreement with 
actual values at the time of overflight. 

We have indicated the difficulties in accurately measuring the soil moisture on the ground 
and thus in establishing the accuracy of any remote sensing techniques. In this report, we 
win assume that the uncertainty in the soil-moisture value for the 0- to 1-cm layer is 0.2 
and the uncertainty for the 0- to 2.5-cm and 0- to 5-cm layers is 0.15. 
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The microwave radiometers used on these flights are listed in table 5. The values for the 
I root-mean-square (rms) temperature sensitivity are those measured over relatively constant 
targets, for example, the Salton Sea or barren desert. Surface temperatures were measured 
diuing some of the flights using a nadir-viewing infrared radiometer (Reference 21) operating 
in the 10- to 12-/ttm atmospheric window. In addition, a scaiming infrared radiometer was 
used to map surface temperature variations for both series. A nadir-viewing camera was used 
to determine the fli^t path of aircraft. This film coverage was then used to determine the 
specific times that the aircraft was over the target fields so that radiometer data could be com- 
pared with ground measurements. 


Table 5 

Microwave Radiometers Used in 1972 and 1973 
Soil-Moisture Flights 


Wavelength 

(cm) 

Frequency 

(GHz) 

Polarization 

Nadir Angle 

Sensitivity 

(K) 

21.0 

1.42 

Nadir 

0° 

2.5 (1 s*) 

6.0 

4.99 

Vertical 

38° 

2.0 (1 s) 

6.0 

4.99 

Horizontal 

38° 

10.0(1 s) 

2.8 

10.69 

Vertical 

38° 

3.0(1 s) 

2.8 

10.69 

Horizontal 

38° 

3.0(1 s) 

1.55 

19.35 

Horizontal 

Scanner 

2.0 (47 ms) 

0.81 

37.00 

Vertical 

38° 

2.4(1 s) 

0.81 

37.00 

Horizontal 

38° 

2.4(1 s) 


*Integration time for observation. 


The calibration stability of the radiometers was checked with the data acquired over water 
and desert targets, in particular the Salton Sea and the desert just east of the agricultural 
area in Imperial Valley. The results for the two years’ flights are presented in table 6. 

The blanks in the table indicate that no useful data were acquired for that radiometer during 
that flight. The values are the averages for between 20 and 40 s of data. The results from 
the Salton Sea observations indicate that in 1972, the 21-cm and 1.55-cm radiometers were 
experiencing gain fluctuations, with those for 21-cm radiometers being more severe. The 
calculated values for a smooth water surface using a salinity of 35 p/t (John D. Hess, private 
communication) are listed for a water temperature of 22°C. Atmospheric effects were in- 
cluded, using atmospheric profiles obtained from rawinsondes launched from Tucson, 
Arizona. The 1973 Salton Sea results are in reasonable agreement with the calculated values. 
The observations over the desert areas show a high degree of consistency, and appear reason- 
able with the exception of 0.8-cm results in 1972 which may be about 1 5 to 20 K too high. 
The variation of microwave Tg over the desert should be less than that for the infrared be- 
cause of the greater thermal sampling depth at the microwave wavelengths. These results 
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Table 6 

Brightness Temperatures Over Water and Desert Targets (in Kelvin) 


Salton Sea Results 


Flight 

Time (LST) 

1.42 

Nadir 

4.99 

Vertical 

4.99 

Horizontal 

10.69 

Vertical 

10.69 

Horizontal 

19.35 

Nadir 

37 

Vertical 

37 

Horizontal 

IR 

1 (3/7/72) 

12:09 

104 


97 


119 

131 

187 

148 

20.7 


12:40 

99 


112 


122 

117 

190 

152 

22.4 

2A (3/11/72) 

11:02 

69 

130 



129 

142 

182 

145 

22.6 


11:32 

40 

126 



118 

141 

187 

150 

25.0 

2B (3/11/72) 

17:09 

52 

126 



118 

108 

196 

155 

25.5 


17:35 

79 

132 



123 

no 

197 

156 

25.6 

1 (2/8/73) 

15:27 

96 

144 


138 

102 

119 

169 

129 

16.5 


15:42 

92 

148 


124 

83 

115 

176 

130 

16.6 

Calculated 


96 

133 

93 

138 

98 

128 

170 

128 


(T„^o = 22°C) 















Desert Results 






1 (3/7/72) 

12:23 

274 


261 


274 

287 

310 

294 

41.3 


12:56 

272 


267 


273 

293 

317 

298 


2A (3/11/72) 

11:15 

262 

283 



269 

295 

304 

289 

45.8 


11:44 

267 

284 



270 

299 

310 

294 

45.7 

2B (3/11/72) 

17:20 

273 

281 



273 

289 

315 

296 

33.8 


17:47 

273 

285 



274 

289 

314 

296 

33.3 

1 (2/8/73) 

15:32 

270 

264 


272 

265 

275 

268 

267 

26.5 


15:47 

270 

266 


270 

262 

273 

270 

269 

27.1 


indicate that the greatest uncertainty in our Tg observations would be for the wettest fields 
which show low Tg ('^ 200 K). The greatest problem occurs for the 1972 21 -cm results, 
where ’s for the Salton Sea were as much as 50 K below the expected result. If the errors 
scale linearly from high Tg , we expect that some of the observed Tg ’s over the wet fields 
may be as much as 20 K below their actual values for the second flight. The errors for the 
other radiometers should be less than 5 K for the wet fields. 


RESULTS 

An example of the radiometer data is given in figure 5. Here the infrared and 21 -cm brief- 
ness temperatures are plotted as a function of distance along the north track over Imperial 
Valley. This fliet path is of particular interest because it includes data over the Salton Sea 
and the uncultivated desert east of the agricultural target area. These data indicate the range 
of brightness temperatures to be observed over such a combination of surfaces, that is, 96 K 
over water, ~ 1 80 to 200 K over the wet test fields, and 280 K over the desert and dry fields. 
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The brightness temperatures of the individual fields were determined by averaging the data 
acquired during the 3-s interval that the aircraft was over each field. The average brightness 
temperature for all the beam positions within the field was determined for the 1 .55-cm 
scanner. Because of the scanning it was possible to determine a 1 .55-cm brightness tempera- 
ture for every field on all of the passes. Therefore, a much larger amount of data was ac- 
quired at this wavelength and will be discussed first. 

1.55 cm-Wavelength Radiometer Results 

Brightness temperatures for the individual fields were compared with ground measurements 
of soil moisture. Figure 6 gives the results for the 1 .55-cm scanning radiometer and the 
average soil moisture in the top centimeter for the three flights in 1972 for light soils (sandy 
loam and loam) and heavy soil (clay loam). The values plotted are the averages of the two 
passes over each field. The results for the two passes agreed within 2 or 3 K for the dry 
fields and 5 or 6 K for the wet fields. The standard deviations were 3 to 4 K for the dry, 
and 8 to 10 K for the wet, reflecting the greater variation in soil moisture expected for a 
wet field. The large amount of scatter in the data for the dry fields is the result of the range 
of surface temperature observed during the different flights. For example, there was a 
15 to 20 K decrease in Tg between the noon and 5 p.m. flights over Imperial Valley on 
March 11, 1972. The range of brightness temperature is the same for both soil types and 
there is a clear linear decrease of brightness temperature with soil moisture. The slope is 
less steep for the heavier soils because of the greater range of soil moistures that is possible 
for them. If the soil moisture is expressed as the percent of FC, this difference can be ac- 
counted for as shown in figure 7. Visually, the scatter in the data is somewhat smaller, and 
quantitatively, the correlation coefficient for these data is slightly greater than for the light 
and heavy soils separately. The results from the 1973 flight, presented in figure 8, support 
this conclusion. The scatter at low values of soil moisture is smaller because the data are from 
one flight only, and the magnitude is reduced by 1 5 to 20 K due to the reduced soil tempera- 
tures for the February flight date. Because of the reduced scatter at low soil moistures, it 
can be seen here that the decrease of Tg did not start at zero soil moisture, but at soil mois- 
tures greater than about 25 percent of FC. The linear-regression results indicated in figure 8 
are for data with soil moistures greater than 25 percent of FC. This lack of variation for 
low values of soil moisture is to be expected from the dielectric constant results (figure 3). 
However, the decrease in Tg occurs at lower values of soil moisture than those at which the 
sharp increase in dielectric constant occurs. This may be a response to the increased moisture 
in the soil just below the 0-to 1-cm layer. When Tg is compared with soil moisture for thicker 
layers of the soil, the flat region extends to higher values of soil moistures, for example in the 0- 
to 2.5-cm layer, the flat region extends to 50 percent of FC or approximately the wilt point of 
the soU. In figure 9, this effect is shown in more detail. Here, Tg values for both the 1972 
and 1973 flights are compared with the excess of the soil moisture in the 0- to 2.5-cm layer 
above WP for the soil -and the result is a linear decrease of Tg with a correlation coefficient 
of 0.866. 
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These results confirm the concept that this radiometer is responding to the soil moisture in 
the top few millimeters of the soil, and that the ability of the soil to maintain a wet surface 
layer depends on the amount of subsurface water present— in this case, moisture above the 
WP in the top 2.5 cm of soil. Expressing the soil moisture as a percent FC for the soil also 
compensates for the effect of soil texture. For the remainder of the report, soil moistures 
will be expressed as the percent of FC for the soil. 

21-cm Wavelength Radiometer Results 

As indicated in figure 5, the 21-cm radiometer data displayed a large response to soil mois- 
ture variation. The T^ results versus the soil moisture in the 0- to 1-cm layer for the two 
years’ flights are presented in figure 10. The fact that the brightness temperature of the 
dry fields for the two years had approximately the same value is partially the result of the 
response of this radiometer to the more stable subsurface temperatures and, perhaps more 
significantly, the result of the instabilities observed in the performance of this radiometer in 
1972. In comparing the 21-cm results with the 1.55-cm results, (figures 7 and 8), we find a 
greater range of Tg observed at 21 cm (about 90 K) than at 1.55 cm (about 70 K) which 
results in a 50-percent greater slope or sensitivity to soil moisture at the longer wavelength. 

Here again there appears to be a region of little variation in Tg to approximately 0.25 FC. 
at which point a more rapid decrease occurs. If a regression analysis is done for the data with 
soil moisture above 0.25 FC, the slope of the curve is increased by 20 percent, but the correla- 
tion coefficient is reduced to 0.840 due to the elimination of the dry field data for which 
the scatter is smaller. The effect of comparing the 21-cm results with the soil moisture in 
thicker soil layers is to extend the flat region to higher soil moisture, for example to 0.5 FC 
or approximately the WP of the soil for the 0- to 2.5-cm layer (figure 11). The comparison 
with the soil moisture in the 0- to 5-cm layer, figure 12, indicates little change from the 
0- to 2.5-cm comparison. The poor correlations observed in figures 1 1 and 12 for soil mois- 
tures above 0.5 FC undoubtedly result from the uncertainty in the soil-moisture determina- 
tions as indicated by the error bars. 

Off-nadir Radiometers 

The results for the vertical channel of the 6-cm radiometer are presented in figure 13. The 
range of Tg (50 to 60 K) is about the same as that observed at 1 .55 cm, and the data from the 
two years agree very well. The decrease of Tg begins for soil moistures greater than ~ 0.25 FC. 
Unfortunately, only a small amount of data was acquired by the horizontal ebanne.1 of this 
radiometer so that information concerning polarization effects was unobtainable. 

Data from the 2.8-cm radiometer were acquired for both polarizations during the 1973 flight 
and are presented in figure 14. There is little (~ 5 K) difference between the polarizations 
, for dry soils and a 1 0-to 1 5-K difference for wet soils. These differences are much smaller 
than those predicted for a smooth surface, that is, approximately 50 K when wet and 30 K 
when dry. This decrease in polarization is due to the surface roughness of the field. Burke 
and Paris (Reference 22) have analyzed this problem and suggested that the effect of rough- 
ness is to make the data appear to have been acquired at a smaller nadir angle, for example 
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at a nadir angle of 30° instead of 49° for the data they were studying. They also suggested 
that the average of the two polarizations should be independent of angle and thus of surface 
roughness. When these concepts are applied to the data in figure 14, the differences indicate 
a nadir angle of between 15° and 20° while the correlation for the average of the polarization 
with soil moisture is not improved over that of either channel taken separately. Also, the 
observed range of this average was not as great as that predicted by Burke and Paris (Reference 
22), which may be due to a roughness effect that cannot be accounted for by a shift of nadir 
angle. 

Results for Vegetated Fields 

As noted earlier in this report, the effect of a vegetation canopy could be approximated by an 
absorbing layer above the soil. The effects of this layer were expected to be strongly depen- 
dent on wavelength, with the longer wavelengths affected least, and the experimental results 
proved this to be the case. In figure 15, values of Tg observed for vegetated fields at wave- 
lengths of 21 and 1.55 cm are presented. The results at 1 .55 cm were, with a few excep- 
tions, independent of soil moisture, while those at 21 cm showed only a slight difference 
from those for bare fields. These results were for fields with covers of small grains or alfalfa 
up to approximately 150- to 20-cm high, that is, a light to moderate cover. In comparing 
Tg values at 21 cm for the vegetated fields (figure 15) with the corresponding bare-field values 
(figure 1 1 ) it is noted that the Tg ’s for the vegetated fields are generally lower than the cor- 
responding bare-field values at the same moisture content. We expect that this is due to the 
cooler soil temperatures and the more uniform moisture profiles for the vegetated fields. The 
latter situation occurs because the transpiring plant draws moisture from its entire root zone 
and not just from the surface soil as is the ease for evaporation from bare soils. As a result, the 
surface layer of the soil does not become as dry as that for a bare field. Thus a 21 -cm radiom- 
eter can sense soil-moisture variations on grasslands or other areas of moderate vegetative 
cover. 

DISCUSSION 

Calculations using the layered model were performed using moisture and temperature pro- 
files observed by the personnel at the U.S. Water Conservation Laboratory at Phoenix 
(Reference 23). The soil moisture and temperature profiles were observed at frequent inter- 
vals after a heavy irrigation on March 2, 1971. The data from the same area at same time of 
year were assumed to be reasonable estimates of the situations occurring during the siir- 
craft overflights. It should be noted that the moisture and temperature profiles had been 
obtained from a smooth field, while the microwave radiometer results were obtained from 
rough-surfaced fields. The moisture profiles used in the calculations are presented in figure 
16. The calculations were performed on predawn (6:30 a.m.) profiles and for the early after- 
noon (1 :30 p.m.) profiles of the same day. The corresponding temperature profiles are pre- 
sented in figure 17. These profiles are probably quite representative of the actual situation 
for the 1972 flights. For the 1973 flights, occurring early in February, the actual tempera- 
tures may be somewhat cooler. 
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The dielectric constants used in the calculation at X = 21 and 1 .55 cm were those pre- 
sented in figure 3. These data were chosen because these soils had textures closest to that 
of the Avondale Qay loam soil at the Phoenix laboratory. The dielectric constants were 
expressed by the Unear regression fits indicated in figure 3. They were for X = 21 cm. 

RP = 2.56 + 030 X SM, SM < 1 1 .5% 

= -930+ 130XSM,SM> 11.5% (9) 

IP =0.06X SM,SM<113% 

= _1 .44 + 0.185 X SM, SM > 1 1 .5% 

and for X = 1 .55 cm 

RP = 2.64 + 0.11 X SM,SM<7.1% 

= -2.1 + 0.78 X SM, 7.1 < 19.0% (10) 

= 735 + 0.25 X SM, SM> 19.0% 

IP = 0.08 + SM,SM< 103% 

= -5.1 + 0.59 X SM, 103 <SM< 19.0% (11) 

= 0.26 + 0.29 X SM, SM > 19.0% 

where RP and IP are the real and imaginary parts of the dielectric constant and SM is the soil 
moisture. The calculated values of T^ , the thermal sampling depth (equation 4), and average 
soil temperature (equation 6) for the 1 8 profiles are presented in table 7. The calculated 
Tg ’s are plotted versus the average soil moisture in the 0- to 1-cm, 0- to 2.5-cm, 0- to 5-cm, 
and 0- to 10-cm layers of each profile and in figures 18 (X = 21 cm) and 19 (X = 1.55 cm). 
These results exhibit the same type of behavior as the aircraft observations, that is, an approxi- 
mately linear decrease of Tg versus the 0- to 1 -cm soil moisture, and an essentially flat region 
at low soU moistures when Tg is plotted versus the soil moisture in the 0- to 2.5-cm, 0- to 
5-cm, and 0- to 10-cm layers. These Tg values are plotted as the X’s in figures 10, 11, and 12, 
containing the aircraft results. The agreement at low and moderate soil moistures is very 
good, while at high soil moistures the calculated values are lower than the observed values. 

This is perhaps a result of the combined effects of the surface rou^ness raising the effective 
emissivity of the soil, and the surface moisture for the wet fields not being uniformly high, 
especially for the furrowed fields. The range of Tg calculated for the dawn profiles is about 
10 K greater than the range for the mid-day profiles. For the moist profiles this is due to the 
rewetting of the surface soil from below (dew formation is unlikely i in the dry desert air). 

For the dry profiles, rewetting is not as great a factor. As a result, the Tg ’s calculated for the 
dawn profiles are always lower 'than those of the previous day.l 

The calculated values of Tg at X = 21 cm were rather insensitive to surface temperature varia- 
tion. While surface temperature changed from 8.5°C on the morning of day 37 to 50°C that 
afternoon, there wasrinly a 10-K change in Tg, and part of that change could be attributed 
to the greater moisture content in the morning. This insensitivity to surface temperature is 
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due to the rather large thermal sampling depth at the 21 -cm wavelength. The depth 
(equation 4) for these profiles ranged from 6 cm for the wettest profiles to about 1 5 cm 
for the day 37 profile. For this depth there is only a 3-K difference in the dawn and mid- 
day profiles. The calculation was also performed for a case with constant 2-percent soil 
moisture, assuming a temperature profile similar to that for day 37. The result was Tg = 273 K 
for the constant moisture profile compared with 279 K for the day 37 profile. The drier 
field had the lower Tg , because the greater thermal sampling depth (34 cm) of its constant 
profile enabled the radiometer to respond to the cooler subsurface temperatures. This effect 
undoubtedly contributed to the scatter in the Tg values observed for the very dry fields 
in figures 10, 11, and 12. 

The calculated Tg ’s at 1.55 cm are plotted as X’s in figure 8. Qualitatively the agreement 
with aircraft results is good for low and moderate moistures. There is the gradual decrease 
of Tg out to about 60 percent of FC and then a sharp decrease. However, the model pre- 
dicts a much lower Tg for wet fields, for example, 195 K at 120 percent of FC compared with 
the observed 235 K. The sharp increase of calculated values of Tg occurs between days 4 
and 5 as the profile changes from being relatively uniform to one with a sharp moisture 
gradient below a dry surface layer (figure 1 6). This dry layer, which causes the sharp increase 
in Tg, forms a surface mulch which reduces the evaporation from the surface. 

The calculated diurnal range of Tg at 1 .55 cm was about 30 K, which is much greater than 
the range at 21 cm. This is because the much shallower thermal sampling depth (~ 1 .5 
to 2 cm) at this wavelength causes a greater sensitivity to the surface soil temperature. The 
aircraft data in figure 8 are for an early February flight with a maximum infrared surface 
temperature of 31°C, while calculations for the dry profile were for late March and early 
April days with surface temperatures of 45 to 50°C. Thus, it appears that the aircraft observa- 
tions are approximately 15 to 20 K higher than calculated values would be for February pro- 
files. Similarly, the calculated values for the dry March profile are about 20 K below the 
values observed for dry fields in the March 1972 flights. This increase in effective emissivity 
is probably due to surface roughness effects which have not been incorporated in the model. 
Volume scattering does not appear to be important for the emission from soils. 

It is expected that volume scattering would cause a lowering of the brightness temperature 
below that which would be predicted by assuming a continuous dielectric. This is the op- 
posite of what is observed. The good quantitative agreement between the calculated and 
observed values at 21 cm indicate that the surface roughness effects are not as important 
at this wavelength. 

The aircraft data presented here are very encouraging for the possible use of a 21 -cm radio- 
meter system 'as a sod-moisture sensor, and they are further supported by recent satellite 
results. Radiometers operating at 21 and 2.16 cm were onboard Skylab, and one at 1.55 
cm is still operating onboard the Nimbus-5 satellite. A sample of the data from these three 
instruments is presented in figure 20. The Skylab data are from a report by Eagleman and in 
Reference 24. The 21-cm radiometer had a 3-dB footprint of 1 10 km on the ground, 
and thus was responding to the average soil moisture over a large swath. The data 
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presented here are for a pass 300-km long across north central Texas in June 1973. The 
greater sensitivity of the 21 -cm radiometer to soil moisture is demonstrated by a 45-K range 
of Tg compared with a 15-K range at a wavelength of 2.16 cm and essentially no variation 
at 1.55 cm. If there had been recent rains and a considerable amount of bare ground, the 
1.55-cm data would have contained significantly lowered values of Tg (Reference 19). Thus 
it would appear that the 21 -cm radiometer is responding to soil-moisture variations and it is 
probable that it is responding through vegetative cover. It is interesting to note that the high- 
est values of T„ observed at 21 cm were in the same range (270 K to 280 K) as our aircraft 
observations. The lowest values (230 K) were considerably warmer than our aircraft results 
(~ 190 K). This is undoubtedly due to the large footprint of the Skylab radiometer averaging 
over a range of surface conditions. Using the data presented in figure 15 for vegetated fields, 
an estimate of the average soil moisture observed by the Skylab radiometer is also plotted in 
figure 20. The results indicate that soil moistures ranged from less than 20 percent of FC 
at the northwest end of the pass to about 70 percent of FC at the southeast end. Additional 
studies by McFarland (Reference 25) with the Skylab data have shown a hi^ correlation 
between 21 -cm Tg ’s and the antecedent precipitation index for this area in Texas and Okla- 
homa. 

CONCLUSION 

The results presented here have shown that the surface emissivity of a soil is determined by 
the dielectric properties of the surface soil layer a few tenths of a wavelength thick while 
the thermal sampling depths are much greater. The variation of the relationship between 
soil moisture and emission as a result of soil texture can be accounted for by expressing 
soil moisture as a percent of field capacity for the soil. The capability of the 21 -cm radio- 
meter to sense sod-moisture variations through a moderate vegetation canopy and the 
promising Skylab results encourage consideration of a radiometer operating at this wavelength 
for routine soil moisture observations. A 21 -cm radiometer with 10- X 10-m antenna would 
have 10-km spatial resolution from space shuttle altitudes (400 km). Information at this 
resolution would be very important in the large scale moisture budget calculations used in 
agricultural yield estimates, in global circulation models, and in antecedent moisture esti- 
mates for runoff predictions. 

Further aircraft and field studies are necessary to refine the relationship between microwave 
brightness temperature and soil moisture and to quantify the effect of vegetation in more 
detail. The possible use of shorter wavelength radiometers, for example, 10 cm, should also 
be studied because of the improved spatial resolution they would afford. 
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Figure 1. Schematic representation of the sources for the radiation received by a 

radiometer at an altitude (H). 



Figure 2. Real and imaginary parts of the dielectric constant 
of pure water at a temperature of 20®C. 



VOLUMETRIC WATER CONTENT, cm^ cm* 



Figure 3. Real and imaginary parts of the di- 
electric constant for clay loam soils at wave- 
length of 21 cm (Reference 10) and 1.55 cm 
(Reference 5). 



Figure 4. Volumetric soil water content versus soil water pressure 
potential for four different soils (Reference 3). 
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Figure 5. Thermal infrared (10-to 12-jLim) and microwave (21-cm) brightness temperature versus 
distance for a flight path at the north end of Imperial Valley. The 21-cm brightness tempera- 
tures and measured soil moistures are indicated for several of the fields. 
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Figure 7. 1 ,55-cm brightness temperatures versus percent 
of FC in the 0- to 1-cm soil layer for the 1972 flights. 
The error bars indicate the uncertainty of the soil mois- 
ture measurement. 
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Figure 8. 1 .55-cm brightness temperature versus percent of FC in 
the 0- to 1-cm soil layer for the 1973 flight. The linear regression 
resuits are for the data with soil moistures greater than 25 per- 
cent of FC. 



Figure 9. 1 .55-cm brightness temperatures versus the soil moisture above the wilt 

, point for the soil. 



PERCENT OF FIELD CAPACITY IN 0- to 1-cm LAYER 

Figure 10. 21 -cm brightness temperatures versus percent 
of field capacity in the 0- to 1-cm soii layer. 
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Figure 1 1 . 21 -cm brightness temperatures versus percent of 
field capacity in the 0- to 2.5-cm soil layer. 
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Figure 12. 21 -cm brightness temperatures versus percent 
of field capacity in the 0- to 5-cm soil layer. 
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Figure 14. Brightness temperature for the off -nadir 2.8-cm radiometer versus percent of field 

capacity in the 0- to 1-cm soil layer. I 
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Figure 16. Soil moisture profiles observed at the 
U.S. Water Conservation Laboratory. The curves 
are labeled by the number of days after irrigation. 
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Figure 18. Calculated 21 -cm brightness temperatures for the 
profiles presented in figures 16 and 17 plotted versus the aver- 
age soil moisture in the 0- to 1-, to 2.5-, 0- to 5-, and 0- to 10-cm 


CALCULATED MICROWAVE BRIGHTNESS TEMPERATURES, KELVIN 



Figure 19. Calculated 1.55-cm brightness temperatures for the pro 
files presented in figures 16 and 17. 
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